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The discriminant-analysis method has been appUed to optimize the exotic-beam charge recognition in a projectile 
fragmentation experiment. The experiment was carried out at the GSI using the fragment separator (FRS) to produce 
and select the relativistic secondary beams, and the ALADIN setup to measure their fragmentation products following 
collisions with Sn target nuclei. The beams of neutron poor isotopes around '^*La and ^"'^Sn were selected to study the 
isospin dependence of the limiting temperature of heavy nuclei by comparing with results for stable ^^''Sn projectiles. 
A dedicated detector to measure the projectile charge upstream of the reaction target was not used, and alternative 
methods had to be developed. The presented method, based on the multivariate discriminant analysis, allowed to 
increase the efficacy of charge recognition up to about 90%, which was about 20% more than achieved with the simple 
scalar methods. 
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1. Introduction 
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The ALADIN experiment S254, conducted in 
2003 at the SIS heavy-ion synchrotron at GSI 
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Darmstadt, was designed to study isotopic effects in 
projectile fragmentation at relativistic energies. Be- 
sides stable ^^'*Sn beams, neutron-poor secondary 
Sn and La beams were used in order to extend the 
range of isotopic compositions beyond that avail- 
able with stable beams alone. The secondary beams 
were produced at the fragment separator FRS [1] 
by the fragmentation of primary -'^'^^Nd projectiles 
with energies of about 900 MeV per nucleon in a 
thick beryllium target. The FRS was set to select 
^^'^ha and, in a second part of the experiment, also 
^"^Sn secondary projectiles of 600 MeV per nucleon 
which were then delivered to the experiment. 

In order to reach the necessary beam intensity 
of about 1000 particles/s with the smallest possi- 
ble mass-to-charge ratio A/Z, it was found neces- 
sary to reduce the thickness of the degrader at the 
dispersive focus (S2) of the FRS (see Fig. 1) that is 
normally used for the selection of isotopically pure 
secondary beams. As a result, the beams arriving at 
the reaction target contained a distribution of neigh- 
bouring isotopes together with the nominal ^^"^La or 
-'^"^Sn. The A/Z of each projectile is known from the 
measured position at the dispersive focus (S2) and 
from its flight time over the more than 80 m beam 
line between the plastic scintillator at S8 (FRS exit) 
and the start detector at the entrance of the AL- 
ADIN setup (Figs. 1 and 2). The atomic number Z, 
however, could only be measured for non-interacting 
projectiles with the MUSIC IV tracking detector 
placed downstream of the ALADIN dipole magnet. 




Fig. 1. Schematic FRS beam line setup. The primary 
beam from SIS enters the production target area (TA) 
from the left. S1-S8 - focal planes. S2 - dispersive fo- 
cus, 88 - secondary beam exit to Cave B and C (after 
|http: / /www- w2k.gsi.de/frs-setupp . 

In order to minimize the upstream material, only 
two thin plastic detectors were mounted near the 
reaction target (Fig. 2): the start detector, for mea- 
suring the arrival time of each projectile, and the 
position detector, for measuring the projectile po- 
sition in the plane perpendicular to the beam di- 
rection. No detector with the resolution necessary 
to resolve individual Z's in the Z « 50 region was 
placed there. It will be shown in the following that, 
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Fig. 2. S254 setup in Cave B. The secondary beams from the 
FRS enter the experimental area from the left (see Ref. [2]). 

with a discriminant method, the atomic number Z 
of interacting projectiles can, nevertheless, be recon- 
structed with high resolution by combining the am- 
plitude and time information available from all the 
upstream detectors (S2, S8, start and position). 

2. Discriminant analysis 

Discriminant analysis (DA) is a statistical method 
for predicting a class membership for a set of vec- 
tors of observations based on a set of observations 
for which the classes are known (a training set). For 
the training set, DA searches for a transform - a 
projection onto a lower-dimensional vector space - 
such that the ratio of the between-class distance to 
the within-class distance is maximized, thus assur- 
ing maximum discrimination. The optimal reference 
frame onto which the projection is made, is found by 
solving the generalized eigenvalue problem for the 
scatter matrices (see e.g. [3-5] and their references) 
defined for the training set. For a training set com- 
posed of K classes, each in number, the total, 
St, and the within-class, Sw, scatter matrices can 
be defined using the notion of covariance matrix: 

St = NCti 

K 

Sw = ^ '^feC'^, 
fc=i 

and the between-class scatter matrix, Sb, can be 
simply expressed as a differencJ^: 

Sb = St — Sw ■ 

Here, N = X^sLi ^^'^ total number of vec- 

tors of observations in the training set and Ct and 
C\y are the covariance matrices, built of the vec- 
tors of observations, for the total set and for each 
individual class, k, respectively. 

^ See e.g. [3,4] for a strict definition 
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The optimal reference frame is found as a sub-set 
of eigenvectors v solving the generalized eigenvalue 
problem: 

Sbv = XStv. (1) 

An important property of DA is its dimensional- 
ity reduction. Since the rank ofSs cannot be greater 
than (K—l) and the rank of a product cannot exceed 
the ranks of the components, the optimal frame, 
in which the classes are best separated, is spanned 
by at most {K — 1) eigenvectors. This means, for 
example, that for a 3-class training set considered 
here (see next section) and 189 dimensional vec- 
tors of observations, the problem gets reduced after 
DA to 2 dimensions, such that the structure of the 
original high-dimensional space is preserved in the 
low-dimensional space. The problem can then easily 
be handled with e.g. 2 dimensional histograms. It 
should be noted that, to avoid the singularity prob- 
lem, the initial dimensionality should not exceed the 
number of vectors of observations [4]. 

DA has a broad range of applications: from its pi- 
oneering apphcation in plant taxonomy [6] , through 
applications e.g. in image recognition, marketing, di- 
agnosis in medicine, to those in nuclear physics [5,7]. 
We found it also very well suited for recognizing the 
charge of secondary beams from the GSI fragment 
separator used as projectiles in the S254 experiment 
[2]. In order to minimize the background from the 
reactions up-stream, the ionization chambers, nor- 
mally used in secondary-beam experiments to deter- 
mine the Z of the beam, were not used in S254. The 
charge of the beam entering Cave B could only be 
measured precisely for the projectiles which did not 
interact with the target, using the MUSIC IV track- 
ing chamber positioned down-stream of the target 
(see Fig. 2) . In order not to overload the data stream 
with these non-interacting-beam events, only a se- 
lected fraction of them has been allowed to trigger 
(after scaling down) and these events were used as 
the training set. 

3. Results and discussion 

From realistic simulations of the beam transport, 
it was found that the most promising observable to 
determine the up-stream charge is the position of the 
beam at the target. During the experiment this was 
measured with the position-sensitive plastic scintil- 
lator detector required for the tracking of projectile 
fragments and positioned just in front of the target 



(see Fig. 2). Optionally, some other simple observ- 
ables could be considered as beam charge sensitive 
as, e.g., the energy loss in the plastic detectors up- 
stream (82, S8, start or position, see Fig. 2) or the 
A/Z variation with the Z of the beam. 

All those beam charge correlated scalar obscrv- 
ablcs, together with the "reference" energy loss in 
the MUSIC chamber and the final multidimensional- 
DA projection, have been tested for their charge dis- 
crimination power. Figure 3 presents distributions 
of these observables for the training set of secondary 
beams of Z = 57 ± 1, for a selected A/Z interval 
[2.183, 2.203] containing the ^^^La isotope, and plot- 
ted in the range of ±4 standard deviations from the 
mean of the total spectrmn. The throe histograms 
in each panel represent three classes of the beam 
charge {Z — 56 - dotted histogram, 57 - solid and 
58 - dashed) and are labeled according to the charge 
identification derived from the energy loss in the 
MUSIC chamber. 

The top panel of Fig. 3 shows a "reference" charge 
distribution of secondary beams as measured by 
the MUSIC chamber for non-interacting projectiles 
from the same run. This observable served to label 
the classes of the training set. It is unavailable for 
projectiles that underwent fragmentation in the 
target but, from the experimental point of view, of 
intc;rc;st for analyses requiring a precise knowledge 
of the atomic number Z besides the known A/Z. 

The three middle panels present the correspond- 
ing "charge distributions" obtained by measuring 
the horizontal position of the beam at the target 
(corresponding to the dispersive coordinate of the 
FRS), the energy loss in the position detector, and 
the high resolution A/Z. 

The bottom panel of the figure shows the same 
three classes of charge plotted as projections on the 
eigenvector corresponding to the largest eigenvalue, 
obtained using the DA technique with 189 dimen- 
sional input vectors of observations (see Fig. 4 for 
a corresponding 2-dimensional original projection). 
The components of these input vectors of observa- 
tions (features) contain the information from all the 
up-stream detectors, i.e. the times and amplitudes 
from the S2 and S8 FRS scintillators [1], the start 
detector and from the position detector in front of 
the target. All together 18 raw signals which were 
combined in addition into pair products to maxi- 
mize the DA efficiency (see e.g. [7]), finally resulting 
in 189 features. 

The figure demonstrates that the high charge res- 
olution obtained with the MUSIC chamber is dif- 
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Fig. 3. (Color online) Charge distributions for Z = 57 it 1 
beams of the training set obtained from the (from top to bot- 
tom) energy loss in the MUSIC chamber, horizontal position 
of the beam at the target, energy loss AEplastic ™ the 
position detector, high resolution A/Z and the DA with 189 
dimensional input vectors of observations projected on the 
eigenvector corresponding to the larger eigenvalue (bottom, 
see Fig. 4). The numbers (top-right) specify the discrimina- 
tion power [7] of the observable. The spectra are plotted in 
the range of ±4 standard deviations of the total. The solid 
{Z = 57), dotted {Z = 56) and dashed {Z = 58) histograms 
represent three classes of the beam charge labeled according 
to the energy loss in the MUSIC chamber. 

ficult to be achieved or surpassed with other sim- 
ple observables. However, combining all the possible 



pieces of information in an optimal way using the 
DA method allows to improve the resolution and to 
come much closer to the reference. 

One way to quantify the resolution of an observ- 
able is to determine its discrimination power. It is 
defined as [7]: 

Sb 
St 



A = 



for 1-dimensional vectors of observations (scalar ob- 
servables) and as an appropriate eigenvalue of (1) 
for the multidimensional case. A S [0, 1] and is equal 
to zero when the classes have the same mean values 
and equal to one when, for all classes, all the vectors 
of observations of a given class are equal but distinct 
from those of the other classes. Discrimination pow- 
ers for the tested observables have been specified in 
the top-right corner of each panel of Fig. 3. 

The discrimination among classes becomes even 
better since in actual discrimination procedures we 
consider at least three classes at a time and thus 
deal with at least two-dimensional projections. Fig. 
4 shows the original 2-dimensional scatter plot for 
the three classes of charge considered here which 
was the base for the one-dimensional projection pre- 
sented in the bottom panel of Fig. 3. The distribu- 
tion results from the projection of the original 189- 
dimensional vectors of observations on the two main 
eigenvectors obtained with the DA technique. 
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Fig. 4. (Color online) Projections of the original 189-dimen- 
sional vectors of observations on the two main eigenvectors 
{vi , V2 ) obtained with the DA method for the 3-class train- 
ing set. The three clouds of symbols represent three classes 
of charge of the example training set: Z = 57 (black circles), 
Z = 56 (blue squares) and Z = 58 (red triangles). 

Having a way to label the classes of the training set 
and the DA tool to define a low-dimensional space 
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in which the classes are best separated, one can now 
project the vectors of observations of unknown class 
membership onto the optimal reference frame and, 
using the distributions of the training set, assign to 
them the most probable (or, if possible, average) 
labels. 

In our application the label (charge) assignment 
has been done in two ways. Since the beam charge 
range is quite narrow, we restricted the individual 
charge identification to 5 classes of charges with Z G 
\Znoin - 2, Znora + 2] (with nominal charge Z„o„i = 
57 (50) for the nominal ^^''La (^°^Sn) beam) plus two 



classes with Z < 



2 and Z > Z„ 



2, thus 



altogether 7 classes. In the first approach this defi- 
nition of classes has been subsequently simplified by 
splitting the 7 classes into 5 groups, each group con- 
taining only 3 classes: Zq with Z G [Z„om—'2, Znom + 
2], Z< with Z < Zq and Z> with Z > Zq. This 
allowed, thanks to the dimensionality reduction of 
DA, to base the charge assignment on the usage of 
2-dimensional histograms, such as the one in Fig. 4, 
defined using the vectors of the training set. 

The other approach was based on the original 
7-class definition which lead to a 6-dimensional 
optimal space after DA. The 6-dimensional distri- 
butions of projections of the training vectors were 
then parametrized assuming their Gaussian shapes 
and these parametrizations were used for the charge 
assignment. The first approach, based on the 2- 
dimensional reference histograms, was found to give 
about 2% better efficacy of charge recognition as 
compared to the method based on the Gaussian 
parametrization. 

Fig. 5 shows the proportion of charges recon- 
structed using the DA method, Zoa, as a function 
of the reference charge measured by the MUSIC 
chamber, Z music - The numbers in boxes represent 
the proportions, in % of the total yield, calculated 
for the total statistics collected. They are indi- 
cated if the contribution exceeds 1%. The sum of 
the diagonal elements provides the measure of the 
overall efficacy of the method, i.e. the proportion of 
correctly recognized charges. It amounts to about 
86% in case of the nominal " ^^''La" beam and to 
about 89% for the "^"^Sn" beam. The first and the 
last bins contain additionally the contributions of 
all the smaller and greater charges, respectively. In 
particular, the figure shows that the light impuri- 
ties formed in reactions on the up-stream material 
(the first column) amount to about 4-6% of the 
incoming projectiles. About 50% of them can be 
properly recognized and eventually rejected dur- 
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Fig. 5. Contributions of charges reconstructed using the 
DA-based method (Zj^a) ^■iid the charges measured by the 
MUSIC chamber {Z music)- The numbers in boxes repre- 
sent the contributions in % of the total yield. The first (last) 
bin contains also the contributions of all the smaller (greater) 
charges. 

ing the analysis. The methods based on the simple 
scalar observables have been found to be about 20% 
less efficacious than the one based on DA. 

The ratio of the mass to atomic numbers, A/Z ^ of 
the projectile has been measured with high precision 
from the time of flight and position at the dispersive 
focal plane S2 of the FRS. The calibration has been 
done using the basic equation for identification: 



Bp 



A 



up 



where, i?p is the projectile rigidity and A, Z and (3 its 
mass and charge numbers and velocity, respectively; 
u and e are the atomic mass and charge units. 
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Fig. 6. Distributions of AjZ for "i24La" (top) and "^O^Sn" 
(bottom) nominal beams. The numbers in the top left corners 
specify the average AjZ of the secondary beams for the 
whole experiment. 

Figure 6 presents the measured distributions of 
AjZ for both secondary beams used in the experi- 
ment. The high precision measurement aUows to vi- 
suahze the "fine" structure of the AjZ groups due 
to its Z dependence. The average {AjZ) for the two 
secondary beams are: 2.188 for the nominal " ^^"^La" 
beam and 2.160 for the "^'^^Sn" beam. The average 
{Z) for the two beams are 56.8 and 49.7, respec- 
tively. The averages take into account variations due 
to different FRS settings for different groups of runs 
and are taken over the whole experiment. 

Having determined the AjZ and Z, the isotopi- 
cally pure beams can be selected also for events with 
reacting projectiles. Figure 7 shows the resulting 
composition of the two secondary beams used in the 
S254 experiment. 

The numbers in the histograms specify the con- 
tribution of a given isotope (in %) and are indicated 
only if the contribution is greater than 1%. Gating 
on a given atomic and mass number allows to se- 
lect a 90% pure beam and to also keep away from 
the light impurities from reactions on the up-stream 
material. 

Summarizing, we have demonstrated a high ef- 
ficacy of the secondary-beam charge recognition 
method based on the discriminant analysis tech- 
nique in an experiment without a dedicated Z- 
detector. The method has improved the quality of 
the S254 experimental data by allowing up to 90% 
pure isotopic selection of the incoming projectiles. 

The authors would like to thank P. Desesquelles 
for stimulating introduction to DA and valuable 
comments. This work has been supported by the 
Polish Ministry of Science and Higher Education 
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Fig. 7. Mass versus atomic number of the two secondary 
beams used in the S254 experiment, as obtained from the Z 
reconstruction using the DA method and the high precision 
AjZ measurement. The numbers in the boxes represent the 
contribution of a given isotope (in %) to the secondary beam. 

under Contracts No. 1 P03B 105 28 (2005 - 2006) 
and N202 160 32/4308 (2007-2009). 
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